Dark Matter Results from 225 Live Days of XENON100 Data 
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We report on a search for particle dark matter with the XENON100 experiment, operated at 
the Laboratori Nazionali del Gran Sasso (LNGS) for 13 months during 2011 and 2012. XENON100 
features an ultra-low electromagnetic background of (5.3±0.6) x 10 -3 events (kg day keV cc ) _1 in the 
energy region of interest. A blind analysis of 224.6 live days x 34 kg exposure has yielded no evidence 
for dark matter interactions. The two candidate events observed in the pre-defined nuclear recoil 
energy range of 6.6-30.5 keV nr are consistent with the background expectation of (1.0 ± 0.2) events. 
A Profile Likelihood analysis using a 6.6-43.3 keV nr energy range sets the most stringent limit on the 
spin-independent elastic WIMP-nucleon scattering cross section for WIMP masses above 8GeV/c 2 , 
with a minimum of 2 x 10 -45 cm 2 at 55 GeV/c 2 and 90% confidence level. 
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Astronomical and cosmological observations indicate 
that a large amount of the energy content of the Uni- 
verse is made of dark matter [1] . Particle candidates un- 
der the generic name of Weakly Interacting Massive Par- 
ticles (WIMPs) [2] arise naturally in many theories be- 
yond the Standard Model of particle physics, such as su- 
persymmetry, universal extra dimensions, or little Higgs 
models. The search for these particles continues with a 
variety of experimental approaches [3]. In direct detec- 
tion experiments, one attempts to observe the nuclear 
recoils (NRs) produced by WIMP scattering off nucle- 
ons [4]. The expected signal features a recoil spectrum 
which falls exponentially with energy and extends to a 
few tens of keV only. The expected low event rate re- 
quires large detectors built from radio-pure materials and 



that are capable of identifying and rejecting backgrounds 
from various sources. 

The XENON100 experiment, described in detail in [5], 
uses liquid xenon (LXe) as both WIMP target and de- 
tection medium, with simultaneous measurement of the 
ionization and scintillation signals produced by particle 
interactions in the active volume. The detector is a cylin- 
drical two-phase (gas/liquid) time projection chamber 
(TPC) enclosing a LXe target mass of 62 kg. An ad- 
ditional 99 kg of the same high-purity LXe, optically sep- 
arated from the target volume, is instrumented as an ac- 
tive scintillator veto. The TPC and the veto are mounted 
in a double-walled stainless-steel cryostat, enclosed by a 
passive shield made from OFHC copper, polyethylene, 
lead, and water/polyethylene. The shield is permanently 



2 



purged with boil-off N2 gas in order to suppress radon 
backgrounds. The LXe is kept at the operating tempera- 
ture of about — 91°C by a pulse tube refrigerator (PTR) 
mounted outside the shield. For the run leading to this 
new result the PTR has been in continuous operation for 
a total of ~20 months. This is the first demonstration, 
to our knowledge, of a LXe detector operated over such 
a long period of time. 

The key feature of the XENON100 TPC is its ability to 
reconstruct the energy and three-dimensional coordinates 
on an event-by-event basis. This enables background re- 
duction by fiducial volume optimization, exploiting the 
self-shielding of LXe. An energy deposition in the TPC 
produces both ionization electrons and scintillation pho- 
tons. The electrons, moved from the interaction point 
by a drift field of 530V/cm, are extracted from the liq- 
uid and accelerated in the gas by a ^12kV/cm field, 
producing proportional scintillation light. The amplified 
charge signal (S2) and the direct scintillation signal (SI) 
are both detected by two arrays of 1" -square Hamamatsu 
R8520-AL photomultipliers (PMTs), selected for low ra- 
dioactivity and high quantum efficiency [5] . One array is 
immersed in the LXe below the cathode of the TPC for 
optimal light collection, and one is placed in the xenon 
gas above the amplification gap. 

The z-position of a particle interaction in the TPC is 
reconstructed, with a precision of 0.3 mm (ler), from the 
time difference between the SI and S2 signals and the 
known electron drift velocity. The localized distribution 
of the S2 signal over the top PMTs is used to obtain 
the (x, j/)-coordinates using a Neural Network algorithm 
with an uncertainty < 3 mm (ler) [5, 6]. The spatial 
reconstruction also allows for the rejection of multiple- 
scatter events, such as from neutrons, since WIMPs are 
expected to interact only once. Double-scatter events can 
be separated when their vertices differ by Az > 3 mm. 
Finally, the ratio S2/S1 is different for NRs (WIMPs, 
neutrons) and electronic recoils (ER; /3/7-background) 
and is also used for background discrimination. 

XENON100 is installed at the Laboratori Nazionali 
del Gran Sasso (LNGS) of INFN, Raly, at an average 
depth of 3600 m water equivalent, where the muon flux 
is suppressed by 10 6 [7]. Due to careful material selec- 
tion [8] and detector design [5] , the total ER background 
of XENON100 is - 5 x 10" 3 eventsx (keVocXkgxday)" 1 
(keVee = keV electron-equivalent [9]) in the dark matter 
energy region, before S2/S1 discrimination. 

Compared to the data reported in [6], the new dark 
matter search is characterized by a considerably larger 
exposure and a substantial reduction of the intrinsic 
background from 85 Kr by cryogenic distillation. The 
nat Kr concentration in Xe has been lowered to (19 ± 
4) ppt, as measured in a Xe gas sample from the detec- 
tor using ultra-sensitive rare gas mass spectrometry com- 
bined with a sophisticated Kr/Xe separation technique. 
This is consistent with the (18 ± 8) ppt derived from 



counting the number of delayed (3-j coincidences asso- 
ciated with the 85 Kr beta decay, assuming a 85 Kr/ nat Kr 
ratio of 2 x 10" 11 [10]. 

The data have been acquired under improved elec- 
tronic noise conditions and the hardware trigger was 
changed to a majority design to allow for a reduced 
S2 trigger threshold of > 99% for S2 signals above 
150 photoelectrons (PE). This has been directly mea- 
sured using the method described in [5, 11] and leads to 
virtually no loss of events in the energy region of interest. 

The non-uniform light collection by the two PMT ar- 
rays and the attenuation of the ionization signal by resid- 
ual impurities over the maximum drift gap of 30 cm lead 
to a position-dependent SI and S2 signal response. The 
signals are corrected using maps derived from calibration 
data. The SI light yield is 3-dimensionally corrected in 
cylindrical coordinates (r, 9, z) in order to optimize the 
response very close to the PMTs. The electron lifetime 
r e [5], used to describe the ionization loss by impuri- 
ties in LXe, was measured regularly with a 137 Cs source 
throughout the data taking period. The value increased 
from 374 /its to 611 /xs, with the average being r e = 514 fis. 
The measured drift time td is used to correct the S2 signal 
size for these losses, and an additional correction in (x, y) 
accounts for variations due to photon collection efficiency 
and small inhomogeneitics in the mesh electrodes. The 
width of the S2 signal is also corrected in (x, y, td) such 
that it is independent of these inhomogeneities. For the 
analysis presented here, the maximum size of the latter 
two corrections is 15% and 3%, respectively. The cor- 
rections, including one due to the imperfect drift field, 
are described in more detail [5] . The considerably larger 
amount of ER and NR calibration data taken during this 
dark matter run (48.0 and 2.7 live days, respectively), al- 
lowed for the improvement of the accuracy of most of 
these corrections at the percent-level. 

As in the previous analysis, it was decided a priori to 
use the Profile Likelihood (PL) statistical inference as 
introduced in [12]. Both the signal and the background- 
only hypothesis are tested. An analysis based on the 
maximum gap method [13] with a pre-defined signal re- 
gion (benchmark region), is used as a cross check. 

The energy scale for NRs is derived from the 
SI signal using the independently measured rel- 
ative scintillation efficiency C c s via the relation 
E nl = (SI/ L y )(l/ £ e g)(S ee / S nr ) (see [14] and references 
therein). The £ e ff parametrization of [6] is used. A hard 
cut in the SI acceptance (see Fig. 1) removes nearly all 
events below the last measured £ e ff datapoint at 3keV nr . 
making the extrapolation to lower E nl irrelevant. The 
factors S ee = 0.58 and 5 nr = 0.95 describe the scintil- 
lation quenching due to the electric field and are taken 
from [15]. L y = (2.28 ± 0.04) PE/keV cc is the updated 
response to 122 keV gamma rays as determined from cal- 
ibration measurements using lines above and below this 
value. The interpolation between these lines is performed 
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using the NEST model for scintillation [16]. 

After verification that the SI pulses are not affected 
by electronic noise, the lower energy threshold used for 
this analysis was set to 3PE, corresponding to 6.6keV nr - 
The PL analysis takes into account the expected WIMP 
energy distribution and would not need an upper energy 
threshold, but for practical purposes, a high threshold of 
43.3keV nr (30 PE) was employed. The benchmark region 
is limited to an upper threshold of 30.5keV nr (20 PE) 
chosen to optimize the signal-to-background ratio. Sig- 
nal (NR) and background (ER) events can be distin- 
guished by their different S2/S1 ratio, where only the 
S2 signal detected by the bottom PMTs, S2&, is used 
since it requires smaller corrections [5], The mean of 
the log 10 (S2fc/Sl) band from ER calibration data is sub- 
tracted in order to remove the energy-dependence of this 
discrimination parameter. 

The dark matter data used for this analysis were ac- 
cumulated over a period of 13 months between February 
28, 2011 and March 31, 2012. Besides 3 interruptions due 
to equipment maintenance, the data were acquired con- 
tinuously. Dark matter data taking was otherwise only 
interrupted by regular calibrations using blue LED light 
(for the PMT response), a 137 Cs source (for monitor- 
ing of the LXe purity), and 60 Co and 232 Th sources (for 
ER background calibration). Overall, the duty cycle of 
XENON100 during this dark matter run was 81%. To 
calibrate the response to NRs, data from an 241 AmBe 
neutron source were taken just before the start and right 
after the end of the run. The two measurements are in 
good agreement. 

Periods with with increased electronic noise or very lo- 
calized light emission in the xy-plane were removed from 
the data, as well as periods in which crucial detector pa- 
rameters such as temperature or pressure fluctuated out- 
side of their normal range. This results in a final dark 
matter dataset of 224.6 live days. In order to avoid analy- 
sis bias, the dark matter data was blinded from 2-100 PE 
in SI by keeping only the upper 90% of the ER band, 
thereby masking more than 90% of the signal region. 

In order to identify valid NR candidate events with 
the highest possible acceptance, several classes of cuts 
and event selections are applied to the data. Their ac- 
ceptance is evaluated on NR calibration data, with the 
exception of quality cuts which might have a time de- 
pendence due to changing detector conditions. These 
are tested on the non-blind part of the science data or 
on the ER calibration data. The first class of cuts are 
basic data quality cuts which remove events that show 
either unidentified peaks or an excessive level of elec- 
tronic noise or light (as from a very high energy event 
or a high voltage discharge). Since only single-scatter 
events are expected from WIMP interactions, the second 
class of cuts identifies such events using the number of SI 
and S2 peaks in the waveform and the information from 
the active veto. Conditions on the size of the S2 and the 



requirement that at least two PMTs must observe an SI 
peak ensure that only data above the threshold and well 
above the noise level are considered. Finally, it is veri- 
fied that several quantities associated with the event are 
consistent, e.g., that the width of the S2 signal, affected 
by electron diffusion in the LXe, is consistent with the 
^-position derived from the time difference between the 

51 and the S2. 

The acceptance of most of the cuts is given versus the 
size of the measured SI signal, used to infer the NR en- 
ergy scale. This is subject to large Poisson fluctuations 
due to the low number of quanta involved. Only the 
acceptance of the S2 threshold condition, S2>150PE, is 
given vs. the SI signal before Poisson-smearing since the 

52 signal fluctuates independently from the SI after the 
initial energy deposition. This quantity is not accessible 
through a direct measurement and has been numerically 
deconvoluted from measured data [11]. It leads to an 
acceptance which depends on the recoil energy spectrum 
once converted into the observed SI signal (see Fig. 2 
of [6]). The cuts and their acceptance determination are 
identical or similar to those described in [11] with the 
exception of the cut against electronic noise which has 
been improved considerably. The acceptances are shown 
in Fig. 1. 

The fiducial volume used in this analysis contains 34 kg 
of LXe. The volume was determined by maximizing 
the sensitivity of this run given the accessible ER back- 
ground above the blinding cut. The ellipsoidal shape 
was optimized on ER calibration data, also taking into 
account event leakage into the signal region. A bench- 
mark WIMP search region to quantify the background 
expectation and to be used for the maximum gap anal- 
ysis was defined from 6.6-30.5 keV nr (3-20 PE) in energy, 
by an upper 99.75% ER rejection line in the discrimina- 

Sl |PE] 

5 10 15 20 25 30 




Energy [keVnr] 

FIG. 1: Combined cut acceptance (solid blue). The S2 thresh- 
old cut S2> 150 PE (dashed red) is independent of possible 
fluctuations in SI and has to be applied to the SI spectrum 
before taking into account the SI resolution. It is conserva- 
tively set to zero below 1 PE. For the cross-check with the 
maximum gap method [13], a hard discrimination cut is used. 
Its acceptance to NRs is shown by the dotted green line. The 
lower analysis threshold is 6.6keV n r (3PE) and extends to 
43.3keV m -, whereas the cross-check is restricted to 30.5keV nr 
(dashed and dotted black lines). 
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tion parameter space, and by the lines corresponding to 
S2>150 PE and a lower line at ~97% acceptance from 
neutron calibration data (see lines in Fig. 2, top). 

Both NR and ER interactions contribute to the ex- 
pected background for the WIMP search. The first is de- 
termined from Monte Carlo simulations, using the mea- 
sured intrinsic radioactive contamination of all detector 
and shield materials [8] to calculate the neutron back- 
ground from (a, n) and spontaneous fission reactions, as 
well as from muons, taking into account the muon energy 
and angular dependence at LNGS. The expecation from 
these neutron sources is (0.17^ ' 07 ) events for the given 
exposure and NR acceptance in the benchmark region. 
About 70% of the neutron background is muon-induced. 

ER background events originate from radioactivity of 
the detector components and from /3 and 7 activity of 
intrinsic radioactivity in the LXe target, such as 222 Rn 
and 85 Kr. The latter background is most critical since it 
cannot be reduced by fiducialization. Hence, for the dark 
matter search reported here, a major effort was made to 
reduce the 85 Kr contamination which affected the sensi- 
tivity of the previous search [6] . To estimate the total ER 
background from all sources, the 60 Co and 232 Th calibra- 
tion data is used, with >35 times more statistics in the 
relevant energy range than in the dark matter data. The 
calibration data is scaled to the dark matter exposure by 
normalizing it to the number of events seen above the 
blinding cut in the energy region of interest. The ma- 
jority of ER background events is Gaussian distributed 
in the discrimination parameter space, with a few events 
leaking anomalously into the NR band. These anoma- 
lous events can be due to double scatters with one energy 
deposition inside the TPC and another one in a charge 
insensitive region, such that the prompt SI signal from 
the two scatters is combined with only one charge sig- 
nal S2. Calibration data show that anomalous leakage is 
most likely below ~8PE. The ER background estimate 
including Gaussian and anomalous events is (0.79 ±0.16) 
in the benchmark region, leading to a total background 
expectation of (1.0 ± 0.2) events. 

The background model used in the PL analysis em- 
ploys the same assumptions and input spectra from MC 
and calibration data. Its validity has been confirmed 
prior to unblinding on the high-energy sideband and on 
the vetoed data from 6.6-43.3 keV nr . However, the model 
does not include a population with S2/S1 values below 
the NR signal region extending down to the lowest en- 
ergies with S2<150PE. This population was found only 
after unblinding and might contribute to the background 
at low SI. 

After unblinding, two events were observed in the 
benchmark WIMP search region, see Fig. 2. With en- 
ergies of 7.1keV nr (3.3PE) and 7.8keV nr (3.8 PE) both 
fall into the lowest PE bin used for this analysis. The 
waveforms for both events are of high quality and their 
S2/S lvalue is at the lower edge of the NR band from 



neutron calibration. There are no leakage events below 
3 PE. The PL analysis yields a p-value of > 5% for all 
WIMP masses for the background-only hypothesis indi- 
cating that there is no excess due to a dark matter sig- 
nal. The probability that the expected background in 
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FIG. 2: (Top) Event distribution in the discrimination param- 
eter space log 10 (S2f,/Sl), flattened by subtracting the distri- 
bution's mean, as observed after unblinding using all analysis 
cuts and a 34 kg fiducial volume (black points). A lower 
analysis threshold of 6.6keV nr (NR equivalent energy scale) 
is employed. The PL analysis uses an upper energy threshold 
of 43.3 keVnr (3-30 PE) and the benchmark WIMP search re- 
gion is limited to 30.5keV nr (3-20 PE). The negligible impact 
of the S2>150PE threshold cut is indicated by the dashed- 
dotted blue line and the signal region is restricted by a lower 
border running along the 97% NR quantile. An additional 
hard S2t,/Sl discrimination cut at 99.75% ER rejection de- 
fines the benchmark WIMP search region from above (dotted 
green) but is only used to cross check the PL inference. The 
histogram in red indicates the NR band from the neutron cal- 
ibration. Two events fall into the benchmark region where 
(1.0 ± 0.2) are expected from background. (Bottom) Spatial 
event distribution inside the TPC using a 6.6-43.3 keV nr en- 
ergy window. The 34 kg fiducial volume is indicated by the 
red dashed line. Gray points are above the 99.75% rejection 
line, black dots fall below. 
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FIG. 3: New result on spin-independent WIMP-nucleon scat- 
tering from XENON100: The expected sensitivity of this run 
is shown by the green/yellow band (la /2a) and the result- 
ing exclusion limit (90% CL) in blue. For comparison, other 
experimental results are also shown [19-22], together with 
the regions (la /2a) preferred by supersymmetric (CMSSM) 
models [181. 



the benchmark region fluctuates to 2 events is 26.4% and 
confirms this conclusion. 

A 90% confidence level exclusion limit for spin- 
independent WIMP-nucleon cross sections a x is calcu- 
lated, assuming an isothermal WIMP halo with a lo- 
cal density of p x — 0.3GeV/c 3 , a local circular veloc- 
ity of v = 220km/s, and a Galactic escape velocity of 
v csc = 544km/s [17]. Systematic uncertainties in the en- 
ergy scale as described by the C e s parametrization of [6] 
and in the background expectation are profiled out and 
represented in the limit. Poisson fluctuations in the num- 
ber of PEs dominate the SI energy resolution and are 
also taken into account along with the single PE resolu- 
tion. The expected sensitivity of this dataset in absence 
of any signal is shown by the green/yellow (la /2a) band 
in Fig. 3. The new limit is represented by the thick blue 
line. It excludes a large fraction of previously unexplored 
parameter space, including regions preferred by scans of 
the constrained supersymmetric parameter space [18]. 

The new XENON100 data provide the most strin- 
gent limit for m x > 8GeV/c 2 with a minimum of 
a = 2.0 x 10" 45 cm 2 at m x = 55GeV/c 2 . The max- 
imum gap analysis uses an acceptance-corrected expo- 
sure of 2323.7 kgx days (weighted with the spectrum of a 
100 GeV/c 2 WIMP) and yields a result which agrees with 
the result of Fig. 3 within the known systematic differ- 
ences. The new XENON100 result continues to challenge 
the interpretation of the DAMA [19], CoGeNT [20], and 
CRESST-II [21] results as being due to scalar WIMP- 
nucleon interactions. 
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